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Abstract 

A system of fermions with a short-range interaction proportional to the scattering 
length a is studied at finite density. At any order a n , we evaluate the complete contribu- 
tions to the energy per particle E{kf) arising from combined (multiple) particle-particle 
and hole-hole rescatterings in the medium. This novel result is achieved by simply 
decomposing the particle-hole propagator into the vacuum propagator plus a medium- 
[t ! , insertion and correcting for certain symmetry factors in the (n — l)-th power of the 

in- medium loop. Known results for the low-density expansion up to and including order 
a 4 are accurately reproduced. The emerging series in akf can be summed to all orders 
in the form of a double-integral over an arctangent function. In that representation the 
unitary limit a — > oo can be taken and one obtains the value £ = 0.5067 for the univer- 
sal Bertsch parameter. We discuss also applications to the equation of state of neutron 
matter at low densities and mention further extensions of the resummation method. 

a , 

PACS: 05.30.Fk, 12.20. Ds, 21.65+f, 25.10.Cn 
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1 Introduction and summary 

04 . Dilute degenerate many-fermion systems with large scattering lengths are of interest e.g. for 
modeling the low-density behavior of nuclear or neutron star matter. Also dilute systems of 
ultracold atoms can nowadays be trapped. Because of the possibility to tune (magnetically) 
atomic interactions through so-called Feshbach resonances, ultracold fermionic gases provide 
an exceptionally valuable tool to explore the (non-perturbative) many-body dynamics at 
weak and strong coupling together with the transition from the superconducting to the Bose- 
Einstein condensed state. Of particular interest in this context is the so-called unitary limit in 
which the two-body interaction has the critical value to support a boundstate at zero energy. 
In that special situation the s-wave scattering length diverges, a — > oo, and the (strongly) 
interacting many-fermion system becomes scale invariant. Its groundstate energy is then 
determined by a single universal number, the so-called Bertsch parameter £, which measures 
the ratio of the energy per particle, E(kf)(°°\ to that of a free Fermi gas, E(kfp ^ = 3fcj/10M. 
Here, kf denotes the Fermi momentum and M stands for the (heavy) fermion mass. 

The calculation of £ is an intrinsically non-perturbative problem which has been ap- 
proached in recent years by numerical quantum Monte-Carlo simulations (in a periodic box). 
As the state of the art, at present a value of £ ~ 0.38 seems to emerge from these calculations 
[1], presumably with still debatable error bars due to finite-size corrections etc. The equation 
of state of neutron matter at low densities has also been studied using quantum Monte-Carlo 
techniques [2, 3]. Due to the very large neutron-neutron scattering length, a nn ~ 19 fm, neu- 
tron matter at low densities, p n = kj/3n 2 < 0.05 fm 3 , is supposed to be a fermionic gas close 
to the unitarity limit. The results of a variety of sophisticated many-body calculations are 
summarized in Figs. 3,4 of ref. [3] and these give indications for a value of £ nn ~ 0.5. 

On the other hand effective field theory methods have been used to rederive and systemati- 
cally improve the low-density expansion for a system of fermions with short-range interactions 



[4, 5, 6]. In particular, effective range corrections and (logarithmic) contributions from three- 
particle scattering have been included. In an unpublished paper, Steele [7] has computed 
several higher order contributions proportional to (akf) 3 and (akf) A . These are generated 
either by multiple particle-particle and hole-hole rescatterings in the medium or they arise 
from particle-hole ring diagrams. Based on some numerical evidence, Steele argued that in the 
limit of large space-time dimensions D the hole-hole contributions would be suppressed. Un- 
der this assumption the resummation of the particle-particle contributions becomes possible 
(in the limit D — >■ oo) in the form of a simple geometrical series: — (2akf /3n)[l + 2akf/n]~ 1 . 
(We are choosing the sign-convention such that a positive scattering length a > corresponds 
to attraction.) The Bertsch parameter following from Steele's approximation is ^ St - ) = 4/9, 
surprisingly close to recent quantum Monte-Carlo results. The validity of Steele's arguments 
concerning the expansion in 1/D has been critically reassessed in the more elaborate work by 
Schafer et al. [8]. There it has been shown that if the strong coupling limit a — > oo is taken 
after the limit D — > oo the universal Bertsch parameter is = 1/2 (see eq.(44) in ref. [8]). 

Moreover, the particle-particle ladders (for D = 4) have been resummed in ref. [8] in the 
form of a phase space integral over a geometrical series, — a[l + akfF pp (s, /t)/7r] _1 , and the 
corresponding Bertsch parameter had the (rather small) value ^ pp ^ ~ 0.24. In the same paper 
an analogous expression for the sum of all hole-hole ladders has been given (see eq.(25) in 
ref. [8]). In that case the geometrical series involves two subtractions (hole- hole ladders start 
to contribute at order a 3 ) and the unitary limit a — > oo does not exist for this term. The same 
feature applies to the sum of all particle-hole ring diagrams, which have been discussed in 
section IV of ref. [8]. Beyond these partial resummation results there exist still large classes of 
diagrams with mixed particle-particle and hole-hole ladders which have not been considered 
in ref. [8] (or elsewhere in the literature). Clearly, these should all be resummed in order 
to see whether they contribute in the limit a — > oo to the universal Bertsch parameter £. 
The first mixed (pp and hh) ladder appears at order a 4 and it has been considered by Steele 
[7]. Hammer et al. [9] have carefully checked Steele's calculation (for D = 4) and they find 
differences for just this term. In particular, the numerical value given for it in eq.(16) of ref. [7] 
(middle term) has to be corrected and multiplied by a factor 2. Actually, Hammer et al. [9] 
find a simpler analytical representation for it in the form of a phase space integral over the 
squared particle-particle bubble times the hole-hole bubble, 3F^ (s, k)F pp (—s, k) (see eq.(7) in 
section 3). Besides this leading term, results for mixed particle-particle and hole-hole ladders 
at higher orders a n ,n > 5 are not known. For large n their classification and combinatorial 
counting is already a non-trivial problem. 

The purpose of the present paper is to close this gap. We will evaluate, at any order a n , 
the complete contributions to the energy per particle E{kf) arising from combined (multiple) 
particle-particle and hole-hole rescatterings in the medium. A key to the solution of this 
problem is a different organization of the many-body calculation from the start. Instead of 
treating (propagating) particles and holes separately, we keep them together and measure 
the difference to the propagation in vacuum by a "medium-insertion". The latter involves 
a delta-function for on-shell kinematics and a step-function which restricts momenta to the 
region inside the Fermi sphere. In that organizational scheme the pertinent in-medium loop 
(or in- medium bubble) is complex- valued. The contribution to the energy per particle E{kf) 
at order a n is therefore not obtained directly from the (n — l)-th power of the in- medium 
loop. However, after reinstalling the symmetry factors l/(j + 1) which belong to diagrams 
with j + 1 double medium-insertions, a real-valued expression is retained for all n. Known 
results for the low-density expansion up to and including order a 4 are accurately reproduced 
in our scheme. The emerging series in a kf can even be summed to all orders in the form of a 
double-integral over an arctangent function. In that explicit representation the unitary limit 
a — > oo can be taken straightforwardly and one finds the value £ = 0.5067 for the universal 
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Figure 1: Left: A closed multi-loop diagram with (at least) two medium-insertions repre- 
senting a contribution to the energy density. Right: After opening at the pair of adjacent 
medium-insertions (symbolized by short double-lines) the planar ladder diagram results. A 
wiggly line symbolizes the contact interaction proportional to the scattering length a. 

Bertsch parameter. As an application we discuss the equation of state of neutron matter at 
low densities and finally we mention further possible extensions of the resummation method. 




2 Preparation: In-medium propagator 

We are interested in the equation of state of a (non-relativistic) many-fermion system with 
a short-range two-body interaction proportional to the scattering length a. In perturbation 
theory the interaction contributions to the energy density are represented by closed multi-loop 
diagrams with a certain number of contact vertices. Their evaluation proceeds via Feynman 
rules which introduce a factor Altai/ M for each interaction vertex and the (non-relativistic) 
particle-hole propagator: 



p -p 2 /2M + ie p -p 2 /2M -%e, 

2n5(p -p 2 /2M)9(k f -\p\), (1) 



Po - p 2 /2M + ie 

for an internal fermion-line carrying energy po and momentum p. The (large) fermion mass is 
denoted by M. The identity in the second line of eq.(l) gives the separation of the in-medium 
propagator G(po,p) into the vacuum propagator and a "medium- insert ion" . Note that by 
construction a medium-insertion puts a particle onto the mass-shell po = p 2 /2M and restricts 
its momentum to the interior of the Fermi sphere \p\ < kf. The Fermi momentum kf is related 
to the density by p = gk^/6n 2 . We are considering here only the case with spin- degeneracy 
factor g = 2. Relevant many-body contributions to the energy density come from diagrams 
with at least two medium-insertions. In case of the closed ladder diagram shown in Fig. 1 
this minimal pair of medium- insertions has to placed on adjacent positions of the double- 
ring, for the following reason. After opening one gets a planar ladder diagram whose energy 
denominators are all equal to differences of fermion kinetic energies. Only for this planar 
topology the resulting factors of M (from the energy denominators) will balance the 1/M 
factors from the interaction vertices such that a finite result remains (at each order a 11 ) in the 
non-relativistic limit. The open ladder diagram shown in Fig. 1 comes in all possible variations 
with further medium-insertions on internal lines. Due to the momentum-independence of the 
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Figure 2: The in-medium loop organized in the number of medium-insertions. The reflected 
partner of the middle diagram with one medium-insertion is not shown. The external momenta 
|_Px j2 | < kf are from the region below the Fermi surface. 



contact interaction all these multi-loop diagrams factorize (successive loops have nothing in 
common) and they can therefore be summed together in the form of a power of the in-medium 
loop. 

3 In-medium loop generated by a contact interaction 

The basic ingredient for calculating the energy density arising from ladder diagrams with a 
short-range contact interaction has been identified as the in-medium loop (or the in-medium 
bubble). As enforced by the two minimal medium- insertions the incoming momenta p\ and 
f>2 are from the region below the Fermi surface jpi^l < kf. It is convenient to introduce their 
half sum P = (pi + f>2)/2 and half difference q = (p± — P2)/2. In our ordering scheme the 
in-medium loop is built up from contributions with zero, one, and two medium-insertions: 
B + E>i + B 2 . The corresponding one- loop diagrams are shown in Fig. 2 and their evaluation 
is exhibited now in detail. 

The contribution Bq with zero medium-insertions is the well-known rescattering bubble in 
vacuum. It is normalized relative to the tree-level interaction and therefore it includes only 
one factor 4nai/M. The evaluation of the loop integral proceeds as follows: 

Altai °T dl r d 3 l i 
° ~ "ir i 2^ J (2tt) 3 (p-2 + p| )/4M + /q _ ( p + f)2/2M + ie 

i 

(Pi + Pf )/4M -Iq-(P- l) 2 /2M + ie 
= 47ra/-— - , =- dl 1+ % . )=0 + iaq , 2 

using residue calculus for the energy integral J^dlo and in the last step the rule dl 1 = 
of dimensional regularization has been applied. In other regularization schemes the occurring 
scale- dependent constant —2fi/ir is absorbed into a -1 to define the renormalized (physical) 
scattering length [7, 8]. The resummation of infinitely many rescatterings in the vacuum in 
the form of a geometrical series leads to the unitarized scattering length approximation: 

/ = a{l + ia\q\ + (ia\q\) 2 + ...}= - 1 = , (3) 

L J a i — i\q\ |g | (cot oo — «J 

which implies the relation tan5 = a\q \ for the s-wave scattering phase shift Sq. 

Next, we turn to the term Bi from diagrams with one medium-insertion. There are 
two equal contributions which differ in their representation merely by the sign of the loop 



Figure 3: Integration regions in momentum space for the real and imaginary part of the 
in- medium loop. 



momentum I and together they read: 

d 3 l 1 



Bi = -4^ / , {B(k f -\P-T\) + 6(k f -\P + r\)}. (4) 

The integral over a shifted Fermi sphere is most suitably performed by using spherical coordi- 



nates with the range 0< |/| < \P\y + \J k 2 — P 2 {1 — y 2 ) oi the magnitude The quantity y 

is a directional cosine which covers the full range — 1 < y < 1 since the displacement vector P 
satisfies the condition \P\ < kf (see also the left part of Fig. 3). Carrying out this procedure 
one obtains the following result for the real part: 



RbB 1 = -^R(8,k), (5) 

7T 



with the logarithmic function: 



R(s, k) = 2 + 1[1 + (s + Kf] In ^ + S + K , +kl + {s- Kf] In i±£^ , ( 6 ) 

written in terms of the two dimensionless variables s = \pi +f>2\/2kf and k = \p% —p2\/2kf. 
Since both external momenta p 12 are from inside the Fermi sphere one has the additional 
constraint s 2 + k 2 < 1. It is worth to note that the function R(s, n) in eq.(6) is equal to the 
sum of the particle-particle bubble and the hole-hole bubble, R(s, k) = F pp (s, k) + F pp (—s, k), 
with the unusual feature that the latter is also taken at momenta below the Fermi surface. 
For comparison the particle-particle bubble reads [7, 8]: 



1 + S + K 1 , 9 9x 1(1 + s) 2 — K 2 \ 

F pp (s, K) = l + s-Kln — 1 + - 1 -s 2 - /? ) In l -\ ^ . 7 

1 + S — K\ 2s 1 — S — K 



Finally, we come to the right diagram in Fig. 2 with two-medium insertions. Obviously, it 
generates a purely imaginary contribution. The total imaginary part of the in-medium loop 
has the following representation: 

r dH 

lm(B + B t + B 2 ) = Ana J —— tt 5{l 2 - f) 



(2k) 

{[l-^y-iP-fD] [i_e(*/-|p + r|) 

+ 0( kf -\p-r\)6(k f -\p+T\)}, (8) 



where we have suitably arranged terms with no, one, and two step-functions #(...). The first 
term in eq.(8) of the form [1 — 0(...)][1 — 9 (...)] makes no contribution to the imaginary part 
since the corresponding phase space is completely Pauli-blocked: 2k 2 < (P — I) 2 + (P + I) 2 = 
2{l 2 — q 2 ) +p 2 +p| < 2k 2 . This equation expresses the obvious fact that on-shell scattering 
of two particles from below the Fermi surface into the region above the Fermi surface is not 
possible due to energy conservation. Thus there remains the imaginary part due to the second 
#(...) 9(...) term in eq.(8). After visualizing the occurring product of 9- and (^-functions one sees 
that the imaginary part has a nice geometrical interpretation: namely as \q\ times that part 
of the solid angle of a (centered) sphere of radius \q\ which lies inside the intersection region of 

— * 

two spheres of radius kf with their centers displaced by 2\P\. The corresponding configuration 
of spheres is sketched in the right part of Fig. 3. After inclusion of the appropriate prefactor 
the result for the imaginary part of the in-medium loop reads: 

Im(B + B 1 + B 2 ) = ^=akfI(s,K), (9) 



with the (non-smooth) function: 



k for < k < 1 - s , 



/(S ' K) = < Ul-s 2 -n 2 ) for l-s<n<VT^, m 

where k lies in the interval < k < y/1 — s 2 . It is interesting to observe that the diagram 
with two medium-insertions alone gives twice the total imaginary part. Putting the real and 
imaginary pieces together the complex-valued in-medium loop reads: 

B + B 1 + B 2 = -^{R(s,K)-inI(s,k)}, (11) 

and if the contribution from the diagram with two medium-insertions is taken out, the imag- 
inary part of that same expression changes sign: 

B + B 1 = -^{R(s,K)+tTTl( S ,K)}. (12) 

Actually, this special property of the in-medium loop turns out to be crucial in order to derive 
the correct expression for the energy per particle E{kf) from powers of the (complex-valued) 
in-medium loop. The derivation of the energy per particle is the topic of the next section. 

4 Energy per particle 

In this section we show how the contributions to the energy per particle E(kf) at any order 
a n can be constructed from the in-medium loop and how the emerging series in akf can be 
summed to all orders. Consider the (open) ladder diagram with n contact interactions. It is 
given by the (n — l)-th power of the in-medium loop times a factor 4nai/M. Closing the two 
open fermion-lines introduces an integration over the allowed phase space |pi )2 | < kf. The 
emerging integrand (R — ml) n ~ l for the energy density at order a n would be complex- valued 
and therefore it cannot yet be the correct one. The deficit of the (naive) iteration method 
at this intermediate stage becomes evident if one draws all diagrams with (repeated) pairs of 
adjacent medium-insertions. The corresponding set of diagrams at fourth order a 4 is shown in 
Fig. 4. These (decorated) diagrams have additional symmetry factors which are not respected 
by the binomial expansion of (R — ml) n ~ l = [(R + ml) + (— 2inl)] n ~ 1 . Furthermore, the 
same symmetry factors correct for the overcounting of certain diagrams as introduced by the 
terms of the binomial series. As a result of this combinatorial analysis one has to reweight 



Figure 4: In- medium diagrams contributing to the energy density at order a 4 . Taking into 
account the proper symmetry factors, their total sum is given by the real-valued expression: 
(R + ml) 3 + (R + ml) 2 (-2ml)(l + 1/2) + (R + ml)(-27:il) 2 + (-2mIf/A = R(R 2 - 
7r 2 / 2 ). Diagrams with a single medium-insertion on the inner or outer fermion-line between 
consecutive interactions are not shown. 



in the binomial series expansion of [(R + iirl) + (— 2ml)] n ~ 1 the j-th power of —2ml coming 
from the diagrams with repeated double medium-insertions with the appropriate symmetry 
factor l/(j + 1). This crucial amendment leads to the following summation formula: 

n— 1 /„ i \ i i 

£(i? + i7r/r 1 ^(-2z 7 rlW . — - = -— -{{R + mlY-{R-mlY}. (13) 



3=0 



Note that the identity ("^/(j + 1) — introduces the binomial coefficients for the 

n-th power of a sum, and in this way one can easily reproduce the result on the right hand 
side of eq.(13). In the diagrammatic representation of the energy density, L-ljO is the number 
of different possibilities to attach j + 1 double medium-insertions on a ring with n segments 
and the additional factor 1/n comes from the n rotations which transform the ring into itself. 

Inspection of the right hand side of eq.(13) shows that the resulting homogeneous poly- 
nomials of degree n — 1 in R and it I are manifestly real for all n. The first five terms read: 
n = 1: 1, n = 2: R, n = 3: R 2 - n 2 I 2 /3, n = 4: R(R 2 - vr 2 / 2 ), n = 5: R A - 2R 2 n 2 I 2 + tt 4 J 4 /5. 
In Fig. 4 and the appended caption it is shown explicitly how the real- valued expression 
R(R 2 — 7T 2 I 2 ) results from the sum of (complex- valued) diagrams once their symmetry factors 
are taken into account. The same graphology gives at second order a 2 : {R+mI) + {—2mI)/2 = 
R, and at third order a 3 : (R + ml) 2 + (R + ml)(-2ml) + (-2mI) 2 /3 = R 2 - 7r 2 / 2 /3. 

At this point we have achieved a representation which includes, at any order a n , the 
complete contributions from all mixed particle-particle and hole-hole ladders. One can even 
go further and sum up the whole series of ladder diagrams to all orders. The pertinent series 
akf(R±i7rI)/7r] n /n can be solved easily by a (complex) logarithm. After inclusion of 
the kinetic energy and the Fock exchange-term (which reduces the Hartree term discussed so 
far by a factor 1 — 1/g = 1/2) the complete expression for the energy per particle reads: 



k ) i 3 _ 48 h. ,fr^ , n ^ n „ ak f ^ k ) 



E(kf) = —fri dss / c^Karctan \ , '— , -}. (14) 



The occurring arctangent function refers to the usual branch with odd parity, arctan(— x) = 
— arctanx, and values in the interval [— 7r/2, it/2]. Other branches of the arctangent function 
are excluded by the weak coupling limit a — > 0, which has to give zero independent of the sign 
of the scattering length a. Note that there occurs a discontinuity (by an amount — it) when 
the denominator passes through zero from positive to negative values. This happens for all 
positive values of akf since the function R(s, re) has a logarithmic singularity at s — 0, re — 1. 
The integral-representation for E(kf) given in eq.(14) has some similarity with the expression 
one obtains from the resummed particle- hole ring diagrams (see eq.(29) in ref. [8]). In that 
case the arctangent function involves two subtractions (particle-hole ring diagrams start to 
contribute at order a 3 ) and the integral extends over a 4-momentum transfer J °° dq Q J °° dqq 2 . 
In the present case the reduction to a double-integral (over the quarter unit-disc) has been 
obtained by employing the following master formula for integrals over the interior of two Fermi 
spheres: 



/ ^Mr F ^ ^ = ^? j ds s " / dKK 7 («, K ^ F ^ K ) > ( 15 ) 

|pi, 2 |<fc/ 71 o 

where s = \p\ + f>2\/2kf and re = \p\ — p2\/2kf. Surprisingly, the imaginary part function 
I{s,k) defined in eq.(10) occurs here as the pertinent weighting function. Integrals over the 
outer region of two Fermi spheres can be reduced in a similar way: 

,6 1 



i vT- 



\pi,2\>k f rfZ 



X 



[s 2 + k 2 - 1) 6(1 + s - re) + 2sk 0( K -1- s) . (16) 



This formula is useful e.g. for evaluating the contributions to the energy per particle which 
arise from multiple hole- hole rescatterings in the medium [7, 8]. 

5 Expansion in powers of akf 

Several orders in the low-density expansion of the energy per particle E(kf) are known [4]. 
The terms stemming from particle-particle and hole-hole ladders can be used as a check of 
our calculation which is organized differently by not treating separately particles and holes. 
We find up to and including fourth order: 

-0.0755733 a 3 k 3 f + 0.0524813 a 4 kj + . . . j . (17) 

As it must be, the linear and quadratic coefficients agree analytically. Since only double- 
integrals are involved the other numerical coefficients can be obtained with high precision. 
One finds again good agreement with existing calculations [4, 7, 9]. The third order coefficient 
is: 0.0861836 - 0.0106103 = 0.0640627 + 0.0115106, where the numbers on the left hand side 
correspond to our separation R 2 — tt 2 I 2 /3 and those on the right side refer to the sum of 
two-fold particle-particle and two-fold hole-hole rescatterings. The same comparison for the 
fourth order coefficient gives: 0.0671902 - 0.0147089 = 0.0383116 - 0.0006851 + 6 • 0.0024758, 
with our separation R 3 — Rn 2 ! 2 versus the sum of triple particle-particle, triple hole-hole, and 
combined particle-particle and hole-hole scatterings. At this point it is important to note that 
a factor of 2 is missing in Steele's numerical result for the latter contribution (see middle term 
in eq.(16) of ref. [7]). This error has been confirmed by Hammer et al. [9] who have carefully 



checked Steele's calculation. It has become evident that the squared imaginary part I(s, k) 2 
of the in-medium loop includes also important many-body correlation effects. Their precise 
mapping into the (traditional) particle-hole counting scheme is not obvious. 

As mentioned in the introduction, Steele [7] has suggested a resummation to all orders in 
form of a simple geometrical series: 

E(M (St) = i(- ^M- (18) 

1 /; 2M15 3n + 6ak f { V ; 

Although the original arguments (via a 1/ D-expansion) have been (partly) invalidated [8] , it 
may serve as a useful reference, in particular since the associated Bertsch parameter £( St ) =4/9 
comes out fairly realistic. We adapt the expansion of the energy per particle in powers of akf. 

^{f+gf(-W}- « w » 

such that all coefficients c n become 1 for Steele's resummation result. Performing the same 
expansion with our complete expression for E(kf) given in eq.(14) we find the following values 
for the first dozen expansion coefficients c n : 

ci = 1, c 2 = ^(11 -2 In 2) = 0.8240319119, 
35 

7T 2 9 

c 3 = — + + 5c 3 = 0.8787170548 , 5c 3 = 2.1365 • 10~ 8 , 

3 12 160 ' 

c 4 = 1.22717534 - — ( 10 - n 2 - 4 In 2 2 + — In 2^ = 0.958529 , 
70 V 3 J 

c 5 = 1.14589, c 6 = 1.37081, 

c 7 = 1.76240, c 8 = 2.19993, 

c 9 = 3.03120, c 10 = 3.74458, 

Cll = 5.85642 , C12 = 5.96732 . (20) 

One observes that the first few coefficients stay below 1 while the higher ones show a tendency 
to increase appreciably with n. The apparent region of "convergence" of the power series in 
akf can therefore roughly be estimated as kf < 1.3/ \a\. Despite various attempts we have 
not succeeded to derive the exact value of c% in terms of (reasonable) mathematical constants. 
The value C3 ~ 7r 2 /12 + 9/160 represents an extremely accurate analytical approximation, 
which is far more precise than all previous determinations of this coefficient [4, 7] . 

For large scattering lengths a the low-density expansion is of very limited validity. In 
that case one has to take the full expression for the energy per particle E(kf) as given by 
eq.(14). The solid line in Fig. 5 shows its ratio to the (free) Fermi gas energy 3/c 2 /10M as a 
function of the dimensionless parameter akf. The behavior of this ratio is demonstrated for 
both signs of the scattering length (a > for attraction and a < for repulsion). Outside the 
region kf\a\ > 6 the curve is almost flat while a sharp peak develops inside with a maximum 
value of about 1.62 at akf ~ —0.9. The dashed curve in Fig. 5 corresponds to Steele's [7] 
resummation via a geometrical series eq.(18). The behavior on the attractive side (a > 0) 
is quite similar while on the repulsive side (a < 0) the artificial pole at akf = —n/2 causes 
essential differences. 

6 Unitary limit 



The unitary limit a — > 00 is of special interest since in this limit the strongly interacting 
many-fermion system becomes scale invariant. The energy per particle is then determined by 




Figure 5: Energy per particle E(kf) divided by the Fermi gas energy 3k 2 /10M as a function 
of the dimensionless parameter akf. An attractive (repulsive) contact interaction corresponds 
to a positive (negative) value of the scattering length a. 

just a constant multiple of the (free) Fermi gas energy: 

E(k ,)(») = (21) 

with £ the so-called Bertsch parameter. Returning to the expression for the energy per 
particle E(kf) given in eq.(14) one sees that the unitary limit a — > oo can be performed 
straightforwardly. The formula for calculating the Bertsch parameter £ reads: 



i 



f=l-— fdss 2 f d/c/c arctan 7 = 0.5067. (22) 

71 J J R(s, K) 

v ' 

The resulting numerical value £ = 0.5067 is to be compared with £( pp ) ~ 0.237 obtained by 
Schafer et al. [8] from the resummation of particle-particle ladders. Actually, ^ pp ^ is defined by 
the principal value integral £^ = 1 — 80 J^ds s 2 {-q 1 ~ s1 dn k I(s, /t)F p ~ 1 (s, k) and the necessity 
to treat the pole-singularity limits the precision in the numerical computation of this number. 
One observes that the additional mixed particle-particle and hole-hole ladders increase the 
Bertsch parameter £ by more than a factor 2. Moreover, the divergence of the subset of 
hole-hole ladders (for a — > oo) encountered in ref. [8] has disappeared. Clearly, the value 
£ = 0.5067 as obtained here via an analytical calculation rooted in perturbation theory is still 
considerably larger than the result £ ~ 0.38 from recent quantum Monte-Carlo simulations 

[!]■ 

In the work by Haussmann et al. [10] a self-consistent, thermodynamically consistent treat- 
ment of the unitary Fermi gas at finite temperatures has been presented. In their non- 
perturbative approach the exact one- and two-particle Green functions serve as an infinite 
set of variational parameters and extensive numerical work enters into the solutions of the 
stationarity constraints and the thermodynamical potentials [10]. The resulting Bertsch pa- 
rameter at zero temperature T = was found to be £ ~ 0.36. This value is remarkably close 
to the experimental determinations by Bartenstein et al. [11], £ = 0.32 ±0.11, and Bourdel et 
al. [12], £ = 0.36 ±0.15. 
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Figure 6: Energy per particle of neutron matter versus the neutron density p n = kj/37r 2 . The 
dash-dotted line stems from the sophisticated many-body calculation of ref.[15]. 

7 Application to neutron matter and outlook 

As an application of our analytical result eq.(14) for the complete resummation of in-medium 
ladder diagrams we consider the equation of state of neutron matter. Due to the very large 
neutron-neutron scattering length a nn = (18.95 ± 0.40) fm [13, 14] neutron matter at low 
densities is supposed to be a Fermi gas close to the unitary limit. Recent quantum Monte Carlo 
simulations [2, 3] based e.g. on the Argonne nucleon-nucleon potential give some indication 
for such a behavior. 

In Fig. 6 we show the energy per particle E(kf) of neutron matter as a function of the 
neutron density p n = kj/37r 2 . The solid line results from our analytical formula eq.(14) 
inserting a = a nn = 18.95 fm for the scattering length and M = M n = 939.57 MeV for the 
fermion mass. The dash-dotted line corresponds to the sophisticated many-body calculation 
by the Urbana group [15], to be considered as representative of realistic neutron matter 
calculations. The dashed line in Fig. 6 reproduces Steele's suggestion eq.(18) in the form of a 
simple geometrical series. One observes good agreement up to rather high neutron densities of 
p n ~ 0.2 fm -3 , where the dimensionless parameter a nn kf reaches values up to a nn kf ~ 34.3. At 
higher neutron densities repulsive effects from three-nucleon forces (which are included in the 
Urbana calculation [15]) start to play a more significant role. The inclusion of the effective 
range r nn = (2.75 ± 0.11) fm for s-wave rm-scattering can also lead to sizeable changes in 
the equation of state, as demonstrated in ref. [16]. Note however, that only the resummed 
particle-particle ladders have been considered in that work. 

Fig. 7 shows the neutron matter equation of state in a different representation by plotting 
the ratio of the energy per particle E(kf) to the (free) Fermi gas energy 3k 2 /10M against the 
dimensionless parameter a nn kf. The dots in this figure reproduce results of various (quan- 
tum Monte- Carlo) calculations of low-density neutron matter and have been taken over from 
Figs. 3,4 in ref. [3] (see original references therein). The solid and dashed curves are the same 
as in the right part of Fig. 5 only continued further out in the parameter akf. 

The complete resummation of ladder diagrams (with a short-range contact interaction 
proportional to the s-wave scattering length a) as achieved in this work suggests numerous 
possible extensions. The s-wave effective range parameter and the p-wave scattering volumes 
should be included via 0(p 2 ) terms in the contact interaction. Asymmetries with respect to 
the spin and/or isospin degrees of freedom can be studied via an appropriate modification of 
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Figure 7: Energy per particle of neutron matter divided by the Fermi gas energy Sk'j/IOM. 
The dots representing various (quantum Monte-Carlo) calculations are taken from ref.[3]. 



the medium-insertion. The generalization of the resummation method to finite temperatures 
would be an equally interesting project. Work along these lines in progress. 

Acknowledgements 

I thank J.W. Holt, A. Schwenk and W. Weise for many useful discussions. I thank H.W. 
Hammer for communicating to me results from own unpublished work which have been very 
valuable. This work is partially supported by the DFG Excellence Cluster "Origin and Struc- 
ture of the Universe" . 

References 

[1] M. MacNeil Forbes, S. Gandolfi, A. Gezerlis, "Resonantly interacting fermions in a box", 
cond-mat/101 1.2197. 

[2] J. Carlson, J. Morales, V.R. Pandharipande, D.G. Ravenhall, Phys. Rev. C68, 025802 
(2003). 

[3] A. Gezerlis, J. Carlson, Phys. Rev. C81, 025803 (2005). 

[4] H.W. Hammer, R.J. Furnstahl, Nucl. Phys. A678, 277 (2000). 

[5] R.J. Furnstahl, J.V. Steele, N. Tirfessa, Nucl. Phys. A671, 396 (2000). 

[6] R.J. Furnstahl, H.W. Hammer, Ann. Phys. (NY) 302, 206 (2002). 

[7] J.V. Steele, "Effective field theory power counting at finite density", nucl-th/00 10066. 

[8] T. Schafer, C.W. Kao, S.R. Cotanch, Nucl. Phys. A762, 82 (2005). 

[9] H.W. Hammer et al., unpublished; and private communications. 
[10] R. Haussmann, W. Rantner, S. Cerrito, W. Zwerger, Phys. Rev. A74, 023610 (2007). 
[11] M. Bartenstein et al., Phys. Rev. Lett. 92, 120401 (2004). 
[12] T. Bourdel et al., Phys. Rev. Lett. 93, 050401 (2004). 
[13] D.E. Gonzales Trotter et al., Phys. Rev. C73, 034001 (2006). 
[14] Q. Chen et al., Phys. Rev. C77, 054002 (2008). 

[15] A. Akmal, V.R. Pandharipande, D.G. Ravenhall, Phys. Rev. C58, 1804 (1998). 
[16] A. Schwenk, C.J. Pethick, Phys. Rev. Lett. 95, 160401 (2005). 



